//- 



NASA Contractor Report 187130 

V 

/ 


On Representation of Mechanical 
Behavior and Stereological 
Measures of Microstructure 


E.T. Onat and S.I. Wright 

Yale University 

New Haven, Connecticut 


June 1991 


Prepared for 

Lewis Research Center 

Under Grant NAG3-834 


NASA 

National Aeronautics and 
Space Administration 

i r? . ii R K C T N T A T I 0 N : ; r N9 i-L 

y C • ! •• * • ! i r A L : M A v I I w A N . > T •_ R L H LOG 1 1 r L 

AV-J*’ • ' r O : M'~ yC T 0:0- Mn^l G'i, ? ur t 
v v I • J-i v. ) M ’ C-tL RJK Unci j3 

9 j j 9 i 0 L 0 b i \> 




ON REPRESENTATION OF MECHANICAL BEHAVIOR 

AND 

STEREOLOGICAL MEASURES OF MICROSTRUCTURE 


E.T. Onat and S.I. Wright 
Yale University 

Department of Mechanical Engineering 
New Haven, Connecticut 06520 


Summary 

Macroscopic homogeneity of a heterogeneous body is defined from various points of 
view. The applicability of the principle of Delesse to a single macroscopically 
homogeneous body is discussed. It is then seen that a function derived from a 
consideration of the area fraction of a phase can serve as a measure of clustering of 
particles of that phase. 


1. Introduction 

A fundamental task of the engineer is to predict the behavior of a structure composed of a 
given material under extreme operating conditions. The structure may be a piece of material 
undergoing a forming operation or a turbine disk or a part of a space vehicle. 

With this task in mind, phenomenological and microscopic studies are conducted on the 
mechanical behavior of the material of interest. The material will, as a rule, exhibit heterogeneities at 
various scales of magnification. For the sake of simplicity consider here the case where the material 
can be taken to be a heterogeneous continuum. Let this continuum be in a state of equilibrium at 
temperature 0o- 

To conduct an isothermal phenomenological (thought) experiment we may extract a spherical 
element of radius R from the material and apply to its boundary a time dependent homogeneous 
deformation F(x) while keeping the boundary temperature at 0o. Here x > 0 denotes the time and F 
is a linear invertible map on R 3 with F(0) = I. The deformation applied to the boundary will cause the 
development of stresses within the element and these stresses will vary from point to point because 
of the assumed heterogeneity of the material. Let a(x) denote the average of this field over the 
element at time x (Note that a(0) = 0). 

The relationship <t> that exists, for a given element of the material, between the histories of 
applied deformation F(x) and the observed average stress a(x) is of great interest to the engineer 
and material scientist. 



The applied deformation will be accompanied by changes that take place in the internal structure of 
the element. Indeed, it is these changes that impose the relationship 0 on the mechanical behavior. Hence, 
the need for microscopic studies of internal structure at various stages of a phenomenological experiment. 

It must be emphasized that the relationship 0 will depend, because of the assumed heterogeneity of 
the material, not only on the constituents and internal structure of the material but also on the radius R and 
the location xc of the center of the element that was extracted from the material: 

O = 0(x c , R ) 

Suppose that the relationship 0(xc, R) ceases to depend on xc as R increases and tends to the limit 0o- In 
other words, if, for every e > 0, there exists a minimal radius R E such that 

R> R E implies p(0(xc, R), 0o) < £ 

where p is a function that measures the “distance” between two constitutive relationships 1 , then we say 
that the material is macroscopically homogeneous from the point of view of mechanical behavior. 

If for a desirable £ the minimal radius R E is much smaller than some of the characteristic lengths 
associated with the structure, loading, etc., one can, for purposes of structural analysis, replace the actual 
heterogeneous material with one which is homogeneous and exhibits the constitutive relationship 0o. On 
the other hand, in some materials such as manufactured composites the size of the elements of the hard 
phase is usually much larger than, say, the grain size so that R E may, in some cases, exceed characteristic 
lengths associated with the structure and loading. In such a case the classical continuum analysis of the 
structure needs to be supplemented by appropriate modifications. 

The next task is to develop constitutive equations that reproduce, with various levels of accuracy, 
the “homogenized” mechanical behavior 0o. The engineer will then choose from this set constitutive 
equations that are appropriate to the problem at hand and will use these equations, together with the 
conservation laws of mechanics and appropriate numerical techniques to obtain the fields of average stress, 
strain, temperature, etc. over the structure. 


1 For a definition of distance between behaviors see Geary and Onat (1974). 
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Most of the recent work on the (approximate) representation of mechanical behavior <X>o can 
be put in the following “canonical” form (cf. Geary and Onat (1974)) based on the notion of state and 
orientation: 

a(t) =/(S(0) /: Z — > V 2 , I c R n 

( 1 ) 

^ = g(S,D) + T a S (t) g: E x V 2 ^ R n 

where suitable initial conditions and invariance requirements imposed on the functions of/ and g are 
omitted for simplicity. 

It should be pointed out immediately that (1) is appropriate for isothermal deformations of 
non-aging solids and, as noted above, is valid only for material elements that are sufficiently large. 

In (1). <7(0 denotes the average Cauchy stress carried by the material element at time t and 
T| denotes the space of symmetric second rank tensors. D and Q are, respectively, the tensors of 
rate of deformation and rotation at that time defined as follows: 

F(0F ! (0 = D + £2 (2) 

where F is the homogeneous deformation applied to the boundary of the element and the dot denotes 
derivative with respect to time. 

S(0 stands for the N parameters that measure those aspects of the present microstructure 
within a sufficiently large material element ( R > R E ) that are relevant to an (approximate) 
description of the future mechanical behavior of the element. We refer to S as the internal state and 
orientation of the material ; although, we recognize that S cannot describe all the details of the 
microstructure if N is not extremely large. 

The first equation above simply states that the current stress follows from the knowledge of 
the current state and orientation S of the material. Most materials of interest suffer changes in their 
internal structure as they deform. The second equation of (1) is a statement as to how these 
changes take place. Note that S, the rate of change of the state and orientation of the material, 
depends only on the present S and on the current values of the rates of deformation and rotation, D 
and £2. Moreover, the dependence of S on £2, as is implied by the second equation in (1), is of a very 
special kind. An important part of the law of evolution for S comes from a consideration of the role 
played by superimposed rigid body rotations. It is then seen (Geary and Onat (1974)) that the N 
parameters that measure the state and orientation of the material element arrange themselves to 
form irreducible tensors: 
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s = (qi, q2» • ■ ■ . q n ). 


Irreducible tensors serve as measures of symmetry, and tensorial state variables with rank 
higher than zero make their appearance when the material is anisotropic or when it develops 
anisotropy during the course of deformation. 

Clearly the engineer would like to know the answers to the following questions : 

(i) For a given macroscopically homogeneous material, what is the minimum radius 
R £ which renders the mechanical behavior of a spherical element independent, to 
within an error e, of the location and size of the element ? 

(ii) What aspects of the present microstructure strongly influence the future 
mechanical behavior and how can these aspects be incorporated into the 
representation (1)? 

In the present paper, we begin a study of these questions for the simple case of a two phase 
material where each phase is composed of a homogeneous continuum. (It will be appreciated that 
such a material may be a satisfactory model for certain composites, granular media, and damaged 
continua). 

Clearly, the mechanical behavior of such a material will depend not only on the behavior of its 
constituents but also on the volume fraction and on the geometrical arrangement of these 
constituents. For instance, consider the classical case where one of the phases is stiff and the other 
is ductile. If the microstructure is such that one can find planes that do not intersect the hard phase 
then the contribution of the hard phase to the stiffness of the composite could be negligibly small. On 
the other hand, in granular media the elements of the hard phase may align themselves to form 
connected substructures during the course of deformation and this would give rise to a great increase 
in the stiffness of the material. 

A purpose of the paper is to investigate the relevance of stereological measures of 
microstructure (volume, area and line fractions) to the questions raised above. Another purpose of 
the paper is to reconsider the celebrated principle of Delesse of the field of Stereology without using 
the notion of an ensemble and the ergodic hypothesis. 

2. Stereological Measures of Microstructure 

We consider the simple case of a planar body that extends to infinity. The body is composed 
of two homogeneous phases which we shall call black and white. We shall assume that the black 
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phase occupies an open set of R 2 and the white phase is assigned to the complement of this set. An 
example of such a body is shown in Fig. 1. 



Figure 1. Two phase body sampled 
by a line segment and a circular area. 

Two stereological measures of the microstructure of this planar body are the area and length 
fractions, say, of the black phase. These are respectively A\JA and l\Jl where A* is the total area of 
the black phase found in a circle of radius R centered at C, A = n R 2 and /*, is the total length of the 
black segments found on the segment AB of length / (cf. Fig. 1). 

Let C be the mid-point of this segment and 0 the angle the line segment makes with the 
horizontal axis x. Clearly, for a given body, the length fraction is a function of the arguments shown 
below 


|=/i(jcc,0;i) (3) 

where xq is the position vector of C. Similarly, 

^=/2(*c;*) (4) 

The functions fa and fa are easily measurable signatures of the microstructure. 

In a macroscopically homogeneous body the area fraction fa will cease to depend on xq and 
will tend to the constant value ab (the area fraction of the black phase) as R increases. In more 
formal language, we say that a body composed of two phases is macroscopically homogeneous from 
the point of view of its constituents, if there exists, for any x c and e > 0, a minimal R£ > 0 such that 
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R>R £ implies \At /A - < e where ab is a constant. Thus the area fraction of the black phase 

within a circular observation window differs from the area fraction ab of the material by less than an 

error of e, if the window has a radius R E or larger. 

For a given material, R £ will be a function of e: 

Re = r (E) 

We shall see below that two macroscopically homogeneous microstructures that have the 
same area fraction ab may give rise to functions ri(£) and T 2 (e) that differ from each other by an order 
of magnitude, for a given e. 

We point out that the previously introduced R £ is different from R £ in that R £ renders the 
mechanical behavior independent of Xc to within an error of e, whereas R E does the same for the area 

fraction. However, we expect that R E will be related to R E in the following sense: consider two 
materials that have the same area fraction ab. If n(£) » r 2(s) f° r these two materials, then we will 
also have R E <0 »R e < 2 ). Hence, our interest in the function r(e). We shall say more about this 
function in the concluding Section of the paper. 

The line fraction function f\ gives more detailed information on the microstructure of the 
material. For instance the fact that a given microstructure has preferred directions will be discemable 
from the the Q dependence off\. 

The question arises: in a macroscopically homogeneous body will/i also cease to depend on 
its arguments 0 and xq as / increases? A naive interpretation of the principle of Delesse 2 (1847) 
(which states that in a multi-phase material, the volume(area) fraction of a particular phase is equal 
to the area (line) fraction of the phase contained in a planar (line) cut through the volume (area)) 
would suggest that this would be the case and moreover we would have 

lim /i = lim /2 (5) 

/ — > OO R-+ OO 

It will be remembered, however, that the quantities such as volume fraction that appear in the 
principle of Delesse refer to the average of measurements made in the neighborhood of a point over 
an ensemble of samples of the material, whereas in the laboratory and as in the present paper one 
works with a single body and one makes measurements over subareas and segments of the body. 


2 The principle was originally derived for volume and area fractions. Rosiwal (1898) , Thomson 
(1930) and Glagolev (1933) extended it to lower dimensions (cf. Weibel (1980)). 
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Thus, it is not surprising that, as the example considered below shows, the assertion (3) does not 
hold for certain microstructures. However a weaker version of (3), (cf. (9)), is valid for all 
macroscopically homogeneous bodies. 

We shall now discuss, in terms of simple examples, the various issues raised above, such as 
appropriate size of an observation window and the validity and usefulness of the principle of Delesse 
for measurements performed over a single body. 

3. Example: Periodic Checkerboard 

As the first example we consider a two dimensional body which extends to infinity and has 
the periodic structure of a regular checkerboard shown in Fig. 2. The body is generated by the 
repeated application of translations of 2e\ and lej (where e\ and €2 are mutually orthogonal unit 
vectors) to the following cell: 

(Black: when (x,y) e B = (fey): 0 < * and > < 1, and 
ceu ■ ' ( 1 < x and y <2 

(White: when (x,y) e {(x,y): 0 < x and y < 2} - B 

This results in the following checkerboard pattern: 

y 


x 

Figure 2. Checkerboard cut with the lines: 
y = mx, m = 1, 2/3 and 1/3. 

Area Fraction 

We now study, in this two phase material, the dependence of the area fraction At, /A on the 
center xc and the radius R of the measuring circle. 

It is easily seen that as A = n /? 2 increases, At, IA becomes nearly independent of the location 
and the radius of the circle and approaches 1/2. Thus, in the language introduced before, there exists, 
for any e > 0 , a minimal R e ' > 0 such that R > /? E ' implies \A b IA - 1/2 I < e for any x c . 

We are interested in the function r(e) which depicts the e dependence of . Although in the 
present simple case this dependence can be studied analytically, we prefer here a computational 
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approach that can be used for more complicated internal structures. Thus, we calculate 1 A b /A 1/2 I, 
for some (jtc, R ) circle and represent the resulting pair (I A b /A - 1/2 1, R) by a point in Fig. 3a. The 
“cloud” of open circles in this figure represents the results of many such calculations. It is then 
easily seen that the graph of r(e) (with the \A b /A - 1/2 I axis now serving as the e-axis ) will 
coincide with the “right” boundary of this cloud if the graph of the boundary is a monotonically 
decreasing function of e. Otherwise r(e) will exhibit discontinuities and its graph will consist of 

certain monotonically decreasing parts of this boundary. 



\AtlA-\l2\ 

Fig 3.a. Plot of (\A b /A - 1/2 i, R) pairs for periodic checkerboard. 


In this figure the oblique line on the left represents the equation 

R 2 e- 1/2 (7) 
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which, as a simple analysis will show, describes, for square domains of observation, the asymptotic 
behavior of the function r(e) as £ — » 0. The second line is parallel to the first one and is drawn for 
comparison of asymptotic behaviors. We see from Fig. 3a that a measurement with radius R > 5 is 
needed to obtain the area fraction of the black phase to within an error of 0.01. 

Somewhat similar calculations were performed for the random checkerboard considered in 
Section 4 which also has the area fraction 1/2. Fig. 3b represents the results of these calculations. 
The two oblique lines of Fig. 3a are added to this figure for purposes of comparison. We see that a 

much larger radius, namely R > 50, is now needed to estimate the area fraction to within an error of 

0 . 01 . 



10 ' 10 ' 10‘ 3 10" 1 

\AbJA- 1/2| 

Fig. 3b. Plot of (\Ab /A — 1/2 I, R) pairs for random checkerboard. 
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Line Fraction 

Next we consider the line fraction l\Jl of the black phase and study the 9 and l dependence of this 
ratio for the midpoint C fixed at the comer of a black square. As was noted above a naive interpretation of 
the principle of Delesse would lead one to believe that the ratio IjJl should approach the area fraction 1/2 as 

/ increases. However, this is not the case: for instance for the 45° cut 

lyi = 1 for all / (cf. Fig. 2). Similar results hold for certain other cuts. Indeed, if we let m be the slope of 
the line segment ( m = tanf?), we find that for 0 < m < 1 

1 * when m = ^ and p and q are odd and 
P‘3 1 have no common divisor other than one Q 

(°) 

otherwise 

We see from (8) that, in a regular checkerboard, (3) fails to hold for an infinite, but countable 
number of directions. But, it can be shown that for all macroscopically homogeneous bodies (3) can be 
replaced by the weaker relation: 


lim = 


/-> ■ 


In 

1 

2 


f 

\fi(xc,0 ;/)- ab I d0 =0 


Jo 


(9) 


where ab is the area fraction of the black phase. 

We continue to study other aspects of the dependence of l\Jl on 6 and / with the aim of shedding 
further light on the principle of Delesse. 

We observe from Fig. 2 that when m is rational the linear cut has a periodic structure. The period 
L p is given as follows: 

( Vp 2 + <? 2 when when p and q are both even or both odd ^ 

2 Vp 2 + q 2 otherwise 

However, when m is irrational the cut has a quasi-periodic structure, i.e. the black segments 
arrange themselves only quasi-periodically over the cut. It is easily seen that this observation can be 
generalized to higher dimensions. In fact, de Bruijn's (1981) elegant construction of the celebrated quasi- 
periodic Penrose tilings is based on a three dimensional cut of a five dimensional periodic checkerboard 
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Figure 4. l^Jl v. / for fixed 0 and *c- 

Fig. 4. shows the dependence of l\Jl on / for fixed slopes. It will be noticed that lb/ l converges 
to a fixed value as / increases but that the rate of convergence depends on the period L p . 

Let e measure the deviation of l\Jl from 1/2: 

e(m, 0 = |/*//-l/2| (11) 

It is clear that e will depend on the slope m and the length / of the segment of measurement and that 
0 < e(m, /) ^ 1/2 by definition. 

To summarize the above findings, we calculated the function e(m, /) at 200,000 evenly spaced 
points ( m , /)• The results of these calculations are shown in Fig. 5 where the value of e at a given 
point is measured by means of a (logarithmic) scale of grayness, a “larger” value of e corresponding 
to a darker shading of gray. It will be noticed that the horizontal dark lines in this figure correspond 
to some of the exceptional values of m listed in the first line of (8). 


li 



m 



Figure 5. Gray scale plot of e(m, /). (Am=0.001 and A/=2.5) 


Fig. 6 is derived from Fig. 5 and exhibits the m dependence of e for three fixed values of /: viz. 
/ = 10, 100, and 1000. The “graphs” for l = 10 and 100, are composed of 1000 discrete points that 
represent the e values calculated at 1000 evenly spaced points of the m-axis. When / = 1000, the 
function e{m, l = constant ) exhibits such sharp peaks at the points listed in (8) that the discrete 
points of evaluation used in the previous figures do not yield a satisfactory “graph” of this function. 
The e values calculated at some of these special points were added to the figure for a better 
description of the m dependence of e. 
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Figure 6. /=constant sections of e(m, l). 

We are now ready to return to the Principle of Delesse. We raise the question: what gage 
length / should one use in the study of this material in order to obtain a reliable estimate of the area 
fraction from one line fraction measurement? It is clear from (8) and from Fig. 6 that when l = 1000 a 
random choice of the slope is likely to result in a very good estimate of the area fraction. 

In order to be more precise about this statement, we introduce the function Pr [e > /} 

which measures the probability of the error e exceeding a value £ > 0 in an observation performed 
with a line segment of length / and random slope m. For each l and for some £, the value of this 
function was determined from Fig. 6 as follows: first we find the total number of points on the m-axis 
amongst the evenly spaced 1000 points, at which the value of e exceeds ^ and then normalize this 
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number by dividing it by 1000 (cf. the case / = 10 in Fig. 5. where a hint is provided concerning the 
nature of this process). 

Fig. 7 shows the % dependence of Pr [e > t,, /} for three values of /. We see that the 
probability of committing a fixed error ^ by a bad choice of the slope m decreases prodigiously as the 
length of measuring segment increases. 



Figure 7. Increasing the gage length decreases the probability of error 

All calculations described here were made with line segments which had a common starting 
point. The influence of the starting point on similar calculations was studied, but we shall not report 
on these studies because they provide little additional insight. 

4. Example: Random Checkerboard 

In the above example the micro structure exhibited translational and rotational symmetries 
and the cell that defined it was small and simple. Here we consider a square checkerboard that has a 
random microstructure and contains one million black and white elementary squares of size lxl. 

The elementary squares are assigned randomly but in such away that the area fraction of the black 
phase is again 1/2. A sub-region of this random checkerboard is shown in Figure 8. We note that 
a microstructure that extends to infinity can be created by using copies of this square checkerboard. 
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Area fraction calculations for this structure were performed and are summarized in Fig. 3b. 
As noted before, in this case a much larger observation window ( R > 50) is needed to estimate the 
area fraction of the black phase to within an error of 0.01. 

We next calculated the line fraction 4// of the black phase by placing an end of the measuring 
segment at the corner of the random checkerboard. The gray scale plot of the function e(m , /) which 
results from these calculations and, as before, measures the deviation of 4// from 1/2, is shown in 

Fig. 9. 

m 



The gray scale plot of e{m, l) for the random checkerboard might suggest to a casual observer 
that there may be exceptional values of m at which e(m, l) does not vanish as / increases. However, 
as Fig. 10 shows more clearly, there are no such exceptional values in the case of the random 
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e(m,l = 10) e(m,l = 100) e(m,l= 1000) 


checkerboard. 1 
the random che 
is very slow coi 



e 



Figure 11. e(m, l) averaged over m as a function of / 

Fig. 11, which is inspired by (9), provides yet another illustration of this fact. 

5. Conclusions 

We recall that in a macroscopically homogeneous material the area fraction of the black phase 
within a circular observation window will differ from the area fraction a b of the material by less than 
an error of e, if the window has a radius larger than the minimal radius r(e) (cf. Section 2 ). We find 
that the function r( e) summarizes important facts about the microstructure. 

To illustrate this remark, consider, for definiteness, the case where the white phase is 
connected, the boundary of each black element is smooth and a b < 1/2. It will then be seen that the 
diameter of any black element is greater than 2r(l- a b ). 

As the following examples show, the function r(e) also provides a definition and a measure of 
clustering of the elements of microstructure. Consider two macroscopically homogeneous materials 
that have the same area fraction a b and whose black phases are composed of copies of just one disk. 

As noted before the functions rj(£)and /" 2 (e) for these materials may have widely differing values for 
some range of e. 

Reflection will show that the inequality ri(e) » r 2 (e) would imply that the black elements in 
the first material are more clustered. In this regard we refer the reader to the examples discussed in 
Sections 3 and 4 and to Figs. 3a, 3b and 8. Note that, as Fig. 8 shows, a periodic structure 
constructed with the smallest cell achieves minimum clustering. 

We observed in Section 2 that r(e) also provides, in a somewhat definite sense, an estimate 
of the minimum radius R E which renders the mechanical behavior of a spherical element independent. 
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to within an error e, of the location and size of the element. This is a difficult issue that needs further 
work. 

We discussed, with the help of examples, the 6 and / dependence of the function f\ (x c , 0 ; O 
that measures the line fraction of the black phase. These discussions showed that a line fraction 
measurement could provide a reliable estimate of the area fraction (The principle of Delesse ) if the 
gage length / is chosen appropriately. This choice is related to the function r(e). It would be desirable 

to undertake an analytic study of this relationship. 

Line fraction measurements performed at a fixed point C and with fixed / define a function on 
the unit circle of directions. In the periodic checkerboard of Section 3, this function tends to a limit, as 
/ _» « which has the value ab for all directions except at the ones generated by the list in (8) and its 
extension. At these exceptional points the limit function assumes the values given in (8). The 
“spikes” that the graph of this function exhibits (cf. top graph in Fig.6) are a sign that the 

microstructure has translational symmetries. 

It is worth noting that microscopic observations of internal structure of a material often 
produce functions defined on an object such as a unit circle (as in the above case), unit sphere, etc. 
Moreover, the Fourier representations of such functions give rise to coefficients that are irreducible 
tensors (cf. Onat and Leckie (1988)). Thus the tensorial state variables (qi(0, <l2(0»— > 9n(0). 
which describe the changing internal state and orientation for an approximate representation of 
mechanical behavior, will, in general, be related to these Fourier coefficients that come from 
measurements of microstructure. A major problem of mechanics of materials is to develop general 
methods for establishing these relationships. It is also clear that in certain cases the mechanical 
behavior and, hence, qi(r) may not depend at all on the aspects of microstructure studied in a given 
experiment. 

To illustrate these points, consider the simpler task of determining the effective elastic moduli 
of a macroscopically homogeneous body composed of two homogeneous isotropic elastic materials. 
Suppose also that we know aspects of the microstructure such as the area fraction ab, the function 
r(e) introduced above, and, say, the fact that the line fraction f\ does not depend on Xc and 0 for 
large /. It is known that just the knowledge of the elastic properties of the constituents and the 
area fraction ab would enable one to find bounds to to the effective moduli of the composite. But in 
certain cases these bounds cannot be close to each other. Indeed if ab =1/2 and if one of the phases 
is connected and less stiff than the other one, then the switching of materials between phases 
would not change ab but could create a stiffer composite, indicating that the bounds to the effective 
moduli may not be close to each other in this case. Note, however, that the switching of materials 
will, in general, change the function r(e). We may therefore hope that a calculation of the effective 
moduli which takes into account this micro structural information might produce better 
results.Nevertheless we believe that for still better results one will need to introduce into analysis 
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parameters that measure the presence and the strength of connectivity and its strength in 
constituent phases. 
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